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PREFACE 


The importance of impurity diffusion techniques for constructing semi-conductor 
devices is well known, and a computer model for simulating diffusion techniques is very 
desirable because data useful for design or analysis can be generated much more 
rapidly with greater economy and less tedium than when obtaining the data experimentally. 

Tnis report consists of five sections and five appendices. In Section i, the 
elementary classical models for idealized diffusion conditions are discussed, because 
the ideas developed are subsequently used in the more realistic models discussed in 
Section 2. These more practical models do not generally allow analytic solutions, but 
require some type of numerical analysis. The numerical techniques which are used 
are outlined in Section 3 with more details concerning implementation given in Appen- 
dices I and II. Section 4 gives some of the results which have been obtained with the 
computer programs implementing the numerical techniques with the programs given 
in Appendices III and IV. Section 5 deals with the special problems of impurity-rich 
interlayers forming between an oxide and silicon. Appendix V gives a set of computed 
curves for sheet resistance, junction depth, and oxide thickness for different diffusion 
schedules. 
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1. IDEALIZED DIFFUSION MODELS 

Many of the practical techniques for diffusion, such as used for 
the fabrication of source and drain regions for the MOST, utilize a two- 
step diffusion process. The two-step process involves a predeposition 
step and a drive-in step which are described in this and the following 
sections. Practical diffusion schedules require a complex model descript- 
ion ; however, the more complex practical model utilizes simpler concepts 
which are used to describe idealized diffusion conditions. In this section 
the equations used to describe idealized diffusion conditions are given. 

Diffusion processes are usually described by Fick's first and 
second laws of diffusion. 1 Fick's first law defines a parameter called 
the diffusion constant or the diffusivity. In a one-dimensional model, 
it is described by 

F = - D 9N (1-1) 

9 x 

where F is the fluxdensity of particles, N is the concentration of 

3 

impurity atoms (per cm ), and D is the diffusion coefficient 
9 

(cm / sec). This mathematical statement implies that the flux 
of particles is proportional to the gradient of the concentration. Particles 
move from denser to less dense regions. The negative 
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sign indicates that diffusion occurs from regions of high concentration 
to those of Lower concentration. 

The requirement that particles be conserved gives Fick's second 
law. In one dimension it is 


c>N _ _ dF 
b t b x 


(l-2a) 


Substituting equation (1-1) into (l-2a) gives; 



( 1 - Zb) 


The diffusion coefficient is assumed to be constant and equation 


(1-2) becomes 


= D- 
dt 



(1-3) 


Solutions to equation (1-3) are usually used to describe diffusions. 
The solutions are of two general classes, those in which the impurities 
are externally added to the semiconductor and those in which impurities 
that have previously been deposited in the semiconductor are allowed 
to diffuse farther from the surface. 

The first class of solutions are for diffusions called pre- 
depositions. A large concentration of impurities are diffused slightly 
beneath the surface of a semiconductor wafer. There is a heavy con- 
centration of impurities in the ambient surrounding the wafer, and it 
is assumed that the surface concentration remains at a constant value 


equal to the solubility limit of the impurity in the semiconductor. 

20 .3 

For boron in siLicon the solubility Limit is about 4. 0 x 10 cm 
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For this type of diffusion the boundary conditions and initial conditions 


N(0,t) = N q 
N( oo, t) = 0 
N(x, 0) = 0 x^O 


(1-4) 


where N q is the solubility limit. 


The solution to (1-3) is then 
N(x, t) = N 0 ERFcj^] 


(l- 5a) 


where 


ERFC(y)- 1 -JL J exp( - u 2 ) du 


(1- 5b) 


The second class of solutions are for drive-in diffusions. The 
drive-in is performed after a predeposition in order to increase the 
penetration depth and decrease the surface concentration. In order to 
get a solution for the impurity profile after drive-in, the profile 
after the predeposition is approximated by 


N(x,0) = 


N 0 , x£h 


0, x>h 


( 1 - 6 ) 


where N Q h = Q is the total number of impurities in the semiconductor 


after the predeposition. 


Q = N ERFC(u) du 


(1-7) 


It is also assumed that the surface is inpenetrable so that the 


flux of particles at the surface is zero. 


( 1 - 8 ) 
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The second boundary condition requires that the concentration of 
impurities goes to zero as x approaches infinity. The semiconductor 
crystal is assumed to be infinite in the positive x direction. This 
can be assumed if the diffusion Length, 2 s/Dt, is small compared to 
the crystal thickness. 

The solution to equation (1-3) for the drive-in condition can be 

3 

approximated by equation (1-9). 


N(x, t)- __Q exp 

'/Dt TT 


4Dt 


(1-9) 


P-N junctions are formed by diffusing into a wafer that has a 
doping level that is opposite in polarity to the diffusing atoms. (Any 
process for introducing impurities into a semiconductor to produce 
some electrical property is called "doping". ) The initial doping 
level is called the background doping. By definition, the junction 
is formed where the diffusing impurity concentration equals the 
background doping. The junction depth for the impurity profile in 
eauation (1-9) is 


xj =|4Dt In 


Q 


1/2 


NbVDtTf 


(1-10) 


The junction depth is important in bipolar transistors because 
of its effect on the base width. For the double diffused transistor, 
the base width is the difference between the junction depth of the 
base diffusion and the junction depth of the emitter diffusion. The 
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injection efficiency, base resistance, current gain, and other 
parameters of the transistor depend on the base width. Since the base 
width is usually small., the junction depths must be controlled 
accurately. 

Another important parameter is the sheet resistance defined by 
equation (l-ll).for an "N-type", or donor impurity, diffusion. 


R. 


%L n n dx 

1) 


(l-ll) 


In equation (1-11), R g is the sheet resistance, q is the electronic 
charge, is the mobility of electrons, and n is the number of 
electrons. Practically, n can be taken to be the same as N, the 
impurity concentration, in many cases. The base spreading resist- 
ance in bipolar transistors, parasitic resistances in M-O-S 
transistors, and all parameters which depend on the impurity pro- 
file depend on the sheet resistance. It is also important to the 
process engineer since it can be measured easily and is an aid in 


determining some characteristics of the impurity profile. 
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2. PRACTICAL DIFFUSION MODELS 


In this section some of the practical conditions for diffusion 
processes are considered. One of the first major deviations from 
the idealized diffusion process which must be considered is the 
effect of oxidation during a drive-in diffusion. The oxide is sub- 
sequently selectively etched to expose other areas to be diffused. 

Or to expose areas that are to be electrically connected to an in- 
put or output of a device. The oxide that is not etched away is 
used for insulation, and, in some field effect devices, as part of 
the device itself. 

The rate at which the oxide grows on the surface of the 
silicon depends on the temperature, the amount of oxide present, 
and on whether it is grown in a stream or oxygen ambient. The 
growth rate of the oxide is given by equation (2-1). ^ 


dx 0 

dt 


A 


B 


+ 2 • 


x. 


( 2 - 1 ) 


In this equation x Q is the oxide thickness, and B and A 
are the parameters plotted as functions of temperature in Figure (2-1) 
and (2-2) respectively. The parameter B is called the parabolic 
rate constant, and B/A is called the linear rate constant. 





8 



Figure (2-3).. The Effect of a Growing Oxide on (a) a Boron Impurity- 
Profile (b) a Phosphorus Impurity Profile. 

After a wafer has been through the predeposition process, the 

oxide which is grown during the drive-in ideally masks the wafer so 

that the predeposited impurities cannot diffuse back out through the 

oxide. There will be some impurity concentration in the oxide, 

however, because of the segregation effect at the silicon- silicon 

dioxide interface. The effect of the growing oxide on a boron and 

5 

a phosphorus diffusion are shown in Figure (2-3). A constant 
called the segregation coefficient relates the concentrations on the two 
sides of the interface. The segregation coefficient is defined by 
equation (2-2). 

^Equilibrium concentration in the silicon 
Equilibrium concentration in the oxide 


( 2 - 2 ) 



The segregation coefficient of boron is about 0.3, and for 
phosphorus it is about 10.0. From Figure (2-3) it can be seen that 
many of the boron atoms are lost to the oxide causing a depletion of 
impurities near the surface. For phosphorus a small amount of the 
impurity is lost to oxide, so that, as the oxide grows inward, it 
pushes the impurities in front of it. This in turn causes the con- 
centration near the surface to remain at a higher value than is 
predicted by equation (1-9). 

As the oxide grows, it consumes part of the wafer. The ratio 

of the thickness of the silicon consumed to the total oxide thickness 

is a constant called °C , which accounts for the difference in the 

densities of the silicon and the oxide. The numerical value of 

£ 

is approximately 0.45. 

With a growing oxide on the surface of the wafer, two deviations 
from the previous model for drive-in diffusion are immediately 
apparent. There is no longer a convenient stationary coordinate 
system. The surface is moving with a velocity of 0.45 times the 
growth rate of the oxide. Furthermore, because of the escape of 
impurities through the surface, the impurity gradient at the surface 
is not necessarily zero. 

The surface boundary condition is derived by assuming that the 
oxide grows only at the silicon- silicon dioxide boundary, and by 
requiring that particles crossing the boundary are conserved. It is 
useful to redefine the spatial coordinate so that the moving boundary 
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is the new reference, y=-0. This is done by the linear transformation 


y = x - ocX Q . 


(2-3) 


Then a particle which is stationary in the x-coordinate system 

at a point corresponding to y>0 will have an instantaneous velocity 

due to the movement of the origin which is given by equation (2-4). 
, dX 

£*-=-<— 2- (2-4) 

dt dt ' ' 


Therefore there is an instantaneous flux of particles because of the 


moving boundary which is 


dX„ 


F^ = - oC- 


N(y, t), y*0. 


(2-5) 


dt 


From equation (1-1), the flux due to the impurity gradient is 


F 


2 



( 2 - 6 ) 


The total instantaneous flux of particles in the semiconductor is the 
sum of the terms in equation (2-5) and (2-6). 

With the assumption that all oxide growth takes place at the 

boundary, a particle inside the oxide has an instantaneous velocity 

due to oxide growth that is equal to the magnitude of the growth 

rate of the oxide. The velocity is in the negative y direction so the 

flux is the negative of the growth rate. Using the symbol C for the 

impurity concentration in the oxide, and D Q for the diffusion 

coefficient the total flux in the oxide is 

dX 0 

F 0= -D 0 _i^- C(y, t) . (2-7) 

by dt 
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Conservation of impurities at the boundary requires that the 
flux of particles at y=0 be equal to the flux at y=0t Equating the 
flux in the oxide and the flux in the semiconductor at the boundary 
gives 

BN i B c i dX Q 

D Ty|y=o + ' D o-|^|y = o- =< K - oc )-at- N (°' t ) (2 ‘ 8) 

where C(0,t) is replaced by KN(0, t) on the right-hand side of equation 
(2-8), and K is the reciprocal of the segregation coefficient. 

For impurities which diffuse sLowly in the oxide, equation (2-8) 
reduces to 


D- 


BN 

TT 


y=0 


dX 

(K-°C) — N(0,t). 


(2-9) 


For some impurities, such as boron and phosphorous, the diffusion 
in the oxide will be negligible. 

The diffusion equation is written in terms of the new coordinate 
by returning to Fick's second law, equation (l-2a). By substituting 
the total flux of particles in the wafer into equation (l-2a), the new 
equation is 


BN ^D 2 N , , dX o BN 

"r ■ — xJ ■■■"'■ t OC • — • 

d t 2 d t 8 y 


( 2 - 10 ) 


Equation (2-10) accounts for the oxide growth during the drive-in 
diffusion with the surface boundary condition of equation (2-9). For 
diffusion with an external source of impurities in an oxidizing at- 


mosphere, the boundary condition would not be the same. 
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The diffusing impurities must each have a charge of plus or 
minus one electronic charge depending on whether the impurity is 
N or P-type. There is an electric field associated with these charges, 
since there must be an impurity gradient for diffusion to occur. The 
electric field exerts a force on these charged particles which gives 
them a drift velocity, and the product of this drift velocity and the 
concentration gives another component of flux. 

The electric field can be calculated by assuming that the 
semiconductor is in a thermal equilibrium. There are several 
orders of magnitude difference in the diffusion cooefficients of holes 
or electrons and impurity atoms. Therefore, the holes or electrons, 
depending on the type of semiconductor, tend to diffuse away faster 
than the impurity atoms. The electric field prevents this from 
happening, so that the impurity current and the hole or electron 
current are equal. The impurity current is known to be small, 
since the impurities move on the order of a micron per hour. Then 
in order to approximate the electric field, the hole current for a 
P-type semiconductor is set equal to zero. 

°= - q D P^f fq 'p Ep (2 ' ll) 


From equation (2-11), where Dp andfJL^ are the diffusivity and 
mobility of holes respectively, the electric field is obtained as 
given by equation (2-12). 

E= kT dP 


qp dy 


(2-12) 
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In equation (2-12), D /n has been replaced kT/q from the Einstein 
relation. 

The impurity atoms will have a drift velocity which is given by 
equation (2-13 ). 

v=-fj E (2-13) 

where jj, is the mobility of impurity atoms, and the negative sign 
indicates that negatively charged particles move in a direction 
opposite to the direction of the field. 

The flux of particles due to the drift velocity, assuming all 
impurities are ionized, is the product of the concentration of 
impurities and the drift velocity. 

F=-/TEN (2-14) 

If the electric field, the oxide growth, and diffusion are ac- 
counted for, the total flux of particles is the sum of the flux terms 

given in equations (2-5), (2-6), and (2-14). 

dX Q 

F = - D-^ - a: N - UEN (2-15) 

dy dt ~ 


From Fick's second law and equation (2-5), the diffusion equation 
which accounts for oxide growth and field aided diffusion is 


BN 

St 


= D- 


a n 
Sy 2 


+ °c- 


dX 

o 

IF 


5N 

9y 


+ 


ay 


ijJLEN). 


( 2 - 16 ) 


To calculate the hole concentration, p, it is assumed that at 
the high temperatures used for diffusion all impurities are ionized. 



The relationship between the hole and impurity concentrations is 
defined by the following two equations. 
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p •= N + n (2-17a) 

pn = ni (2 -17b) 

Substitution of the electron concentration, n, from equation 
(2-17a) into equation (2-17b) gives an expression for the hole con- 
centration in terms of the impurity concentration and the intrinsic 


concentration, n^. 

N / N 2 2 

P - —x — + I — 7 — + n ; 


1/2 


(2-18) 


) 


If the number of impurities is much greater than the intrinsic 
concentration, the number of holes is approximately equal to the 
impurity concentration. For this condition, the diffusion equation can 
be rewritten using equations (2-12) and (2-15). 

dN 


M. .D-Jiu.ct- 
-» r 


dX. 


at 


+ D- 


8 2 N 


ay 


dt 


ay 


ay ' 1 


(2-19) 


Equation (2-19) shows that the diffusion coefficient is effectively 
doubled for high impurity concentrations under the influence of a 
built in electric field. For low impurity concentrations, the 
number of holes is approximately equal to the number of intrinsic 
carriers, and the field aided term in the diffusion equation is small. 

The electric field adds a term to the flux used in determining the 
surface boundary condition, equation (2-9). For a predeposition 
this term may be significant, because of the high surface concentration, 



but the boundary condition for predeposition is not calculated from 

the flux of particles. For drive-in diffusion the temperature is 

usually much higher, but the intrinsic concentration is also higher. 

The temperature dependence of the intrinsic concentration is given 

7 


by equation (2-20). 

n. = 3.85 x 10 16 T 3/2 exp 


1.21 


2kT 


( 2 - 20 ) 


19 20 

For typical drive-in temperatures, n- is between 10 and 10 , and 

the surface concentration will be lower than the intrinsic concentra- 
tion. Since the hole concentration, equation (2-18), is almost 
constant; the field, given by equation (2-12), is insignificant. 

g 

Lehovec and Slobodskey have investigated field aided diffusion 
by lumping the diffusion and field terms into a single term with an 
effective diffusion coefficient, D' 1 '. 


D* = D(l + f) 


( 2 - 21 ) 


The function f accounts for the difference in the hole concentration 
and the impurity concentration. Using this relation to account for 
the field, they have shown a region of high impurity concentration 
in which the diffusion coefficient was doubled, and region in which 
the field was not significant. 

9 

Equation (2-10) has been used by Grove to model diffusion with 
a growing oxide for wafers with a constant initial doping, and Kato 
and Nishi^ have used this equation to approximate the solution 


for the impurity profile for the more general case. 
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3. COMPUTATIONAL METHODS, 

The practical diffusion models require the application of 
numerical methods to integrate the diffusion equation. This section 
describes some finite difference methods useful for obtaining solu- 
tions. The simplest finite difference method is called the explicit 
method which is discussed in Appendix I. The method consists of 
the repeated use of equation (3-1). 

n‘ + 1 - RNj * N> { 1 - GA 4. R [ <F> -P) )/pj -2 ] +N‘ + , 

.[GA + R . + R <pj +J -Pj ) / Pj + j ]} ( 3 - 1 ) 

In equation (3-1) all variables are the same as defined in Appendix I. 

The difference approximation becomes more accurate as the 
increments, At and Ay. approach zero; but there are further 
restrictions for a solution to be convergent. For equation (1-3), 
the value of R is restricted as is indicated in equation (3-2). 

R = AT 0 < R < 1/2 (3-2) 

AY 2 

Assuming a similar relation holds for the more general 
diffusion equation, equation (3-1) can be used for predicting impurity 
profiles. Starting with the initial condition, the impurity profile is 
predicted at the next time instant using the present profile. For 


each time 
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step, the vaLues of impurity concentration are calculated at the 
discrete intervals iav for all values of j between, but not including, 
zero and J. These two points must be reserved for the boundary 
conditions. This will require J-2 calculations for each time step. 

If the total diffusion time is divided into I increments, the J-2 
calculation must be made 1-1 times to obtain the final impurity 
profile . 

In calculating the impurity profile by the explicit method, each 
predicted value of concentration depends only on three previous 
values as is shown in equation (3-1). The value of R is restricted 
to certain limits for stable solutions. The implicit method, described 
in Appendix II, overcomes both of these difficulties; but the ease of 
computation is lost. A matrix of difference equations must be 
solved simultaneously, but for this method each predicted value of 
concentration depends on all of the previous values and on the 
predicted values. The implicit method should be more accurate 
than the explicit method because of this increased dependence. 

The solution to the diffusion equation is calculated in a way 
that is similar to the explicit method, except that the entire impurity 
profile is calculated by solving the matrix equation. A simplified 
flow diagram that shows the order in which the calculations proceed 
for a computer solution to the diffusion equation is shown in Figure 
(3-1). A complete listing of the program used for predicting solutions 
to the equation by the implicit method is given in Appendix III. 



Read time and temperature 
schedule 


Calculate diffusion coefficients, 
oxide growth rate parameters, 
normalization factors, etc. 


Determine ambient condition 
and set boundary condition 
accordingly 


Calculate terms for matrix equation 
and solve by Gaussian elimination 


/ specified 
diffusion time 
.been reached/ 


llncrease time 


I Ye s 


Calculate junction position, 
oxide thickness, and sheet 
resistance 


Figure (3-1). A simplified flow diagram for the computer solution 
of the diffusion equation 
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The boundary condition at the surface is met by imposing the 
restriction given by equation (2-9) to calculate the surface concen- 
tration. Writing equation (2-9) in finite difference form and solving 
for the surface concentration gives the equation used to meet the 
surface boundary condition. 


N(0, T) = 


r d X -i 

1+AY (K -oc) — / D 

dt 


(3-3) 


The oxide growth rate is calculated from equation (2-1) according 
to the ambient condition. The ambient condition can be changed 
during a diffusion and the boundary condition will still be met. 

The second boundary condition is that the impurity concentration 
goes to zero as y approaches infinity. For a numerical solution, 
some convenient point must be ctiosen where the impurity concen- 
tration is known. If a point is chosen far enough away from the 
surface, the concentration can be specified to be zero. Since the 
impurity profile between the surface and the junction is usually the 
only part that is of interest, the solution will not. be significantly 
affected in this region if the specified point is slightly in error. 

For the computer program given in Appendix III, the junction 
depth is predicted from equation (l - 1 0 ) , and the impurity concen- 
tration is set to equal zero at three times the junction depth. This 
allows the program to handle diffusions for shallow or for very deep 



20 

junctions; although, solutions for shallow junctions are more accurate 


because the increments will be smaller. 


Although the predeposition can be numerically simulated to get 
the initial condition before drive-in, it was found that the impurity 
profile after the drive-in was not affected significantly by the 
deviation from a complementary error function, equation (l-5b), 
during the predeposition. The effect of field aided diffusion is taken 
into account empirically for the predeposition by adjusting the diffusion 
coefficient and assuming that the impurity concentration still 
follows the complementary error function profile. The profile is 
calculated by taking the integral indicated in equation (i-5a), which 
is integrated numerically by the Gaus s - Laguer re quadrature method. 

For the case where a predeposition is not followed by a drive-in, 
such as the emitter diffusion for a double diffused transistor, the 
effect of the electric field cannot be accounted for empirically. This 
type of diffusion can be handled numerically by a slight alteration 
of the program given in Appendix III. 

The sheet resistance is determined by numerical integration 


of equation (3-4). 



G dY 
0 


(3-4). 


In this equation G is the conductivity which is given as a function 

11 

of impurity concentration by Irvin. 



21 


4. RESULTS 

Both the explicit and implicit methods have been used to solve 
the diffusion equation. The explicit program has been run on the MSU 
UNIVAC 1106 and the MSFC Xerox Sigma V. Checks on computational 
accuracy have been made by comparing with the erfc(x) and gaussian 
tables and by comparing with previous work on oxidation done by 
other methods. Some comparisons with experimental data have also 
been made. These results are reported in this section. 

Solutions by the explicit method were found to be accurate 
enough for many diffusion problems. It is well suited for the gener- 
ation of sheet resistance and junction depth data for diffusions in which 
the ambient condition during a drive-in is constant. The implicit 
method also works well for drive-ins with one or two ambient changes 
if the times in each ambient are about equal. If the time in one 
ambient condition is much greater than another, then the shorter time 
will be divided into a smaller number of time increments, and the 
effect of the short ambient change may not show up properly in the 
solution. 

Changing the ambient during a diffusion creates another problem 
in the stability of the solution. As previously stated, the value 
of R, AT/ AY^ , must be chosen between zero and one half in 
order to obtain a stable solution to equation (1-3) by the explicit 
method. If the numerical solution is to be general, then equation (2-16), 
which accounts for the oxide growth and the electric field, must be 
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solved. For this equation the range in which R must be chosen for a 
stable solution may not be the same. For a predeposition, the elec- 
tric field is significant, but there is no oxide so that an equation 
different from equation (2-16) is solved. In the case of a drive-in 
that starts out in an oxidizing ambient, the concentration near the 
surface drops rapidly. The electric field is negligible, in this case, 
after the impurity concentration drops below the intrinsic concen- 
tration and still another equation is solved. 

A single value for R which gives the smallest error for each 
of these conditions cannot be found for the explicit method. The 
implicit method should give better results in general because the 
dependence of solutions on the value chosen for R is greatly reduced. 
Fewer calculations are required to show the effect of an ambient 
change, but more computation is required for each prediction of the 
impurity profile. If one hundred increments in the Y-coordinate 
are used and R is one fourth, the explicit method requires about 
three seconds of execution time on the UNIVAC 1106. Under the 
same conditions, the implicit method requires about thirty seconds 
of execution time; but by choosing R to be unity, better results are 
obtained and the execution time is comparable to that of the explicit 
method. 

The solutions given by Grove were simulated numerically for an 
oxide growing on a wafer with a constant initial background doping 
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Figure (4-1). Comparison of Predicted Results with the Results 
of Grove. 

level. Since reported values of diffusion parameters vary widely, 
the diffusion coefficient and segregation coefficient given by Grove 
was used. The result of this diffusion is given in Figure (4-1). 

The circles represent the extremeties of expe rimentally measured 


concentration. 


distance from surface (cm x 10 ) 


Figure (4-2). Comparison of Predicted Results with the Results 
of Kato and Nishi 

The diffusion data given by Kato and Nishi was also simulated. 
Figure (4-2) is a comparison of the computer solution and the ex- 
perimental measurements. The impurity profile given is for a short 
predeposition followed by a fifty minute drive-in in an oxygen ambient. 
The diffusion coefficient and segregation coefficient used were the 
values reported by Kato and Nishi. 
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Figure (4-3). The Effect of the Electric Field on a Predeposition. 

Figure (4-3) shows a calculated impurity profile with the electric 
field considered and one in which the field term has been removed 
from the diffusion equation. These profiles are for a predeposition 
with boron at 980 degree centigrade for one hour. The most signifi- 
cant difference in these two profiles is that the profile calculated 
with the field-aided term included shows more impurities in the 
semiconductor. The effect of these added impurities as a function 
of the predeposition temperature is shown in figure (4-4). The 
per cent change in sheet resistance and junction depth were 
obtained by comparing the results of computations of impurity 
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Figure (4-4). The Effect of the Electric Field as a Function of 
Temperature . 

profiles with field aided diffusion and with the field equal to zero. 

For these calculations the surface concentration was assumed to be 
constant at 4. 0 x 10^. 

In some cases the ambient is changed several times during a 
diffusion in order to shape the impurity profile to give some desired 
electrical characteristic. Table (4-L) shows four diffusion schedules 
with several changes of ambient and a comparison of calculated and 
experimental values of sheet resistance and junction depth. 
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TIME (MINUTES) 


(a) (b) (c) 

(d) 


N 2 

14 

14 

8 

8 

°2 

15 

15 

10 

10 

h 2 o 

10 

10 

10 

10 

°2 

5 

5 

5 

5 

n 2 

10 

10 

10 

10 

°2 

60 

0 

0 

60 

N 2 

10 

0 

0 

60 

R s (MEASURED) 

160 

140 

150 

200 

R g (CALCULATED) 

188 

149 

172 

206 

Xj (MEASURED) 

Xj (CALCULATED) 

2. 0 

1.3 

1.2 

2. 0 

2. 08 

1.44 

1.2 

1.98 


Table (4-1 ). Some Predicted and Experimental Results 
for Several Ambient Changes. 

Some results which show the effect that the oxide growth has on 
impurity profiles as time progresses are shown in Figures (4-5) 
and (4-6). Figure (4-5) shows that for phosphorus the oxide rejects 
the impurities and pushes them inward keeping the concentration 
near the surface relatively high. For the boron profile, Figure (4-6), 
the concentration near the surface is depleted. A comparison of the 
two figures shows that the boron profile contains fewer impurities 


than the phosphorus profile. 
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5. PREDEPOSITION DIFFUSIONS 

There are several methods used for accomplishing the pre- 
deposition diffusion: (1) a gaeous impurity source such as diborane 

or phosphine is mixed with carrier gases of nitrogen, oxygen, and 
argon; (2) a liquid source such as POCI3 is used with a carrier 
gas mixture which bubbles through the liquid; (3) a solid source 
such as boron nitride wafers or P2O5 is used with a carrier gas; 

(4) a doped oxide source is deposited on the wafer surface either by 
pyrolytic techniques or by "spin-on" methods using a liquid emulsion 
which is dried and baked. 

The first three methods require that the surface concentration 
be maintained at the solid solubility limit in order to obtain control 
of the impurity profiles. High surface concentrations result in lattice 
strains which are not always annealed out by subsequent drive-in 
diffusions. These lattice strains may produce deleterious effects in 
terms of device performance. Consequently, there has recently been 
a good deal of interest in method (4) which can produce predeposition 
profiles with acceptable control and lower surface concentrations. The 
other methods are also widely used and method (1) was studied in this 
work. 

The diborane predeposition procedure is usually carried out in 
a furnace at about 1,000° C. (980° C . in the furnace at MSFC.) 

The carrier gases are N2 , at 4,000cc/min flow rate, O2, at 120cc/min, 
and argon at 22cc/min with 1 percent, diborane in the argon. The 
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diborane argon mixture is supplied for times ranging from about 15 
minutes to 60 minutes. A layer of boron glass (B 2 O 3 ) is formed 
with a thickness up to about 1, 200 Angstroms. This layer is quickly 
dissolved in a 10:1 HF etch. There is an inter-layer which is a 
solid mixture of boron- silicon- oxygen which is resistant to the HF 
etch. The layer ranges from about 100 to more than 200 Angstroms 
according to ellipsometer measurements (which may not be reliable. ) 

Exactly how much of this recalcitrant layer that remains after 
the 10:1 HF etch is dependent upon the time used for the predeposition. 
The data obtained at this time indicates that for a 15 minute pre- 
deposition, the layer is almost completely etched away in HF. The 
surface sheds water readily, and a subsequent short drive-in diffusion 
in does not produce a dramatic drop in the sheet resistance. For 
a 30 or 60 minute predeposition, the layer which remains after an 
HF etch is about 200 Angstroms thick, water clings to the surface, 
and a short drive-in in Ng (5-10 minutes) produces a drop in the 
sheet resistance by a factor of four. An etch cycle in hot (105° C) 
nitric followed by a 10:1 HF ei.ch is effective in removing the layer. 
Three cycles seemed to be sufficient. A short oxidation (5 minutes) 
at 1150° C followed by a 10:1 HF etch is also effective in removing 
the layer. 

Wafers which were predeposited for 30 and 60 minutes were 
divided into two groups. Half' of the wafers were etched (HNO 3 +HF) 
and half were not. These wafers were then subjected to drive-in 
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diffusions in N 2 at 1150OC. for a total time of 50 minutes. The wafers 
were removed at intervals and the sheet resistance checked. Of course 
there is an ambiguity concerning the effective time because of the 
heating and cooling times; however, the results are shown in Figure 5.1 
Apparently, the boron- rich interlayer serves as an effective unlimited 
source for a period of 50 minutes. Since the chemical composition 
of the interlayer is unknown, any calculation of the total number of 
impurities based on the layer thickness must be considered to be 
purely speculative. An estimate based on the assumption that the 
layer is roughly a 50 percent mixture of B 2 O 3 and Si0 2 gives 10l 6 cm ~ 2 
for the effective Q of the source. This is sufficient to maintain a 
saturated surface concentration for roughly an hour at 1150°C . 


40 



Figure 5.1 Sheet Resistance vs. Time for Drive-in Diffusion in 
Nitrogen at 1150° C with Boron- Rich Surface Layer. 
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If the wafers with the boron-rich interlayer are placed into an 
oxidizing ambient, either dry oxygen or steam, the results are quite 
different. The sheet resistance at first holds steady with the drive- 
in time and then begins to increase. The increase will continue until 
such a time that the impurity gradient at the junction is too low to 
maintain a space charge layer. Then the junction becomes "leaky" 
and the resistance measured is no longer related directly to the 
resistance of the diffused P-layer above the junction. Apparently, the 
oxide begins to grow between the boron-rich layer and the silicon, and, 
no doubt, some transformation takes place within the layer. Con- 
sequently, the layer is masked from the silicon, and probably dilated, 
so that it no longer supplies impurities at a rate which maintains a 
saturated surface concentration. 

The best model which one can construct at this time on the 
basis of the available data is as follows. An effective surface source 
with strength Q g ^ (cm ) is assumed for all predeposition times be- 
tween 30 and 60 minutes. The value chosen for Q is approxi- 
mately 3. 5 X 1016 cm For 15 minutes or shorter predep times, 

Q S 1 =0. In order to incorporate this predep model into the overall 
two-step diffusion model, the following scheme is employed: (1) For 

a drive-in in a N 2 ambient, the flux of impurities at the surface is 
calculated and integrated with respect to time. The surface concen- 
tration is assumed to be at the solid solubility limit until the inte- 
grated flux is equal to Q s j. Then the flux is set to zero. If Q s j 
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is incorrect, then the model eventually becomes invalid. However, 
in most practical cases only a short drive-in in N2 is required, so 
that this is not a severe limitation. (2) For a drive-in in an oxi- 
dizing ambient, the surface flux is assumed to be zero until an 
oxide 800 Angstroms thick is grown. Then the flux is assumed to 
follow the law dependent upon the oxide growth rate, the segregation 
coefficient, etc. The model is to some extent arbitrary, and other 
models may give equivalent results. A more rigorous approach 
would deal with the chemical conversion of the interlayer plus the 
oxide growth and diffusion through the oxide. This is obviously a 
very involved problem requiring micro-chemical analysis of thin 
layers for boron, silicon, and oxygen content. Until such empirical 
data are available, the present model is no more nor less speculative 
than any other. 

The curves given in Appendix V for diborane predepositions 
without removal of the boron- rich layer were obtained using the model 
described in the preceding paragraph. The curves are calculated 
using the "explicit" integration program which is given in Appendix IV. 
The program as given has a set of instructions to generate the input 
data. These instructions may be easily replaced with read instructions 
or other instructions to input the data. The program uses the fol- 
lowing data: (1) RSI, the sheet resistance after the predeposition, 

(2) TE1, Til, the predep temperature and time (degree C., and 
minutes,) (3) TE2, the drive in temperature (deg. C), (4) Tl, T2, T3, 
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the drive-in are obtained by assuming that the field- aided diffusion 
during the predep results in an effective diffusion constant of 2D. 
Then the surface concentration is adjusted to give the correct sheet 
resistance for the predep, RSI. This is done by a subprogram 
called FINDNCKRSl, LA, NO). Practically, this usually results 
in a junction depth after the predep which is too deep and gives a 
bias to the junction depths for drive-in. Sheet resistance values 
seem to agree fairly well. The junction depth error is roughly 
. 2 microns for a 60 minute predep. 
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Figure 5.2 Sheet Resistance and Junction Depth vs. Time for 
®2^6 , Diborane, Predeposition Diffusion. 

(Circled point calculated by including time in 0 ) 
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The sheet resistance after the predep can be calculated fairly 
accurately by taking into account the aiding electric field. Curves 
are shown in Figure (5.2) for sheet resistance versus predep time 
at 980°C. It is questionable whether or not it is practical to cal- 
culate the predep profile using the program including the field aided 
effect each time the drive-in profile is calculated. There are several 
diff iculties which must be resolved, one being the problem of recon- 
ciling the two choices of "regions for the solution." If one chooses 
a reasonable region for solving the drive-in problem, this region is 
inevitably too large for an accurate numerical integration for the 
predep profile. On the other hand, choice of two regions, one for 
predep and one for the drive-in, requires a coordinate transformation 
type of operation to translate data calculated using a finer grid to a 
description with a coarser grid. However, one accomplishes this, 
extra computing time and added program complexity are required. 

At the moment it is not certain that the increase in accuracy war- 
rants this approach, but it is being considered. It seems now that 
the most practical approach is to use the sheet resistance data as 
input data to solve the drive-in problem. 

The "explicit" program can also be used to simulate pre- 
depositions using boron nitride wafers for diffusion source s. In this 
case, the predep sheet resistance is calculated from formulas based 
on Goldsmith's (et.ai.) data. JJ For this option, the control index 
ICON is set equal to zero. 



36 


6. CONCLUSIONS AND RECOMMENDATIONS 

The model and program described in this report is useful for 
calculating junction depths and sheet - resistances for practical two-step 
diffusion processes. In cases where the model has been compared 
with empirical data, the results for junction depths have agreed to 
within less than 10 percent, better than the uncertainty of measurement 
for some individual measurements. Sheet resistance data agrees well 
in some cases. Experimental data shows that the sheet resistance 
after a diborane predep is sometimes erratic. Variations of as 
much as 8ohms from the average of 27 ohms for a 60 minute 
diborane predep have been observed. The program cannot, of course, 
cope with this problem. More controllability of the process is needed 
and this requires more understanding of the physics and chemistry 
of the process. When the experimental values for the predep sheet 
resistance are used to predict the value after drive-in, the predicted 
value agrees within 10 percent of the measured value. 

There are some practical effects known to be missing from 
the model. The non-linearity introduced with very heavy phosphorous 
concentrations and resulting in the "emitter -dip" phenomena are not 
modeled in the program. This is not a severe limitation for dif- 
fusions used for MOST devices. There is also a tendency to get 
away from this type of diffusion because of the deleterious effects 
produced as well as the unpredictable nature of the process. The 
complexity of the conversion of the impurity- rich interlayer upon 
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subsequent oxidation has not been modeled. More experimental 
data will be required before a more accurate model can be constructed. 

One feature of the model which is deficient and will be sub- 
sequently remedied is in the method of calculating the sheet resistance. 
The sheet resistance is obtained by integrating the conductivity from 
the surface to the junction. Some of the data in Appendix V are ob- 
tained using the resistivity of the diffused concentration rather than 
the net concentration. In most cases, this does not matter^ but 
for a light doping it: will. The other problem is that the junction 
becomes leaky when the impurity concentration is light and the gra- 
dient at the junction is low. A calculation for the impurity gradient 
at the junction will be inserted to indicate when the calculations are 
unreliable. Practically such a diffusion is useless, but one would 
wish to know when such results will be obtained. 

Additional studies should be made to apply these programs to 
diffusions from doped oxides and from the boron nitride type source. 
Curves can be obtained for the boron nitride predepositions and the 
program can be used to study the P -well diffusions used in CMOS 
work. This problem is very difficult because the P -region must be 
lightly doped and the junction must be deep. Consequently, it is 
very difficult, to produce a junction which is not leaky or "lost" 
altogether because the diffused region becomes N-type. 

Further studies are needed to improve the model for the 
impurity -rich interlayer. In this case, data are needed to learn 
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what is actually happening physically. Chemical analysis of the layer 
before and after oxidation is needed. The program should be modi- 
fied so that the predep profiles calculated using the field- aided model 
are readily incorporated into the computations for the drive-in dif- 
fusion. Subsequent work will integrate the features of the explicit 
and implicit programs and utilize the field- aided predep model if 
the accuracy improvement is felt to warrant the complication. 

At this point the programs can be used to simulate a variety 
of diffusions and for experimentation. Further attention will be 
given to simplifying the procedures for using the programs. This 
will probably result in several programs, each incorporating the 
basic integrated program mentioned above but differing in the input 
data and output data instructions. 



REFERENCES 


1.. A. M. Smith, Fundamentals of Silicon Integrated Device 

Technology , Vol. 1. , ed. Robert M. Burger and R. P. 
Donovan, Englewood Cliffs: Prentice -Hall, 1967, p. 18 7. 

2. Ibid., p. 189. 

3. Ibid. , p. 195. 

4. A. S. Grove, Physics and Technology of Semiconductor- Dev ices, 

John Wiley and Sons, inc. , New York, 1967, pp. 28-29. 

5. Ibid., p. 70. 

6. A. S. Grove and others, "Redistribution of Acceptor and Donor 

Impurities During Thermal Oxidation of Silicon, " Journal 
of Appl ied Ph ysics , September, 1964, Vol. 35, pp. 2695-2701. 

7. W. R. Runyan, Silicon Sem iconductor - Technology , McGraw-Hill 

Book Co. , New York, 1965, p. 164. 

8. Kurt Lehovec and Alexis Slobodskoy, "Diffusion of Charged 

Particles Under Consideration of the Built-in Field," 
Solid-State Elec tronic s, 1961, Vol. 3, pp. 45-50. 

9. A. S. Grove and others, pp. 2695 -2701. 

10. Taketoshi Kato and Yashio Nishi, "Redistribution of Diffused 

Boron in Silicon by Thermal Oxidation, " Japanese Journal 
of Applied Physics , Vol. 3, pp. 377- 383, 1964. 

11. John C. Irvin, "Resistivity of Bulk Silicon and of Diffused 

Layers in Silicon," The Beil System Technical Journal, 
March, 1962, Vol. XL'i, pp~387--4/0. 

12. M. L. Barry and P. Olof sen, "Advances in Doped Oxides as 

Diffusion Sources, " Solid State Technology, October, 1968, 
Vol. 11, No. 10, pp. 39~42. 

13. N. Goldsmith, J. Olmstead, and J. Scott, Jr., "Boron Nitride as 

A Diffusion Source for Silicon, " RCA Re view , pp. 344 -350, 
June, 1967. 


39 



BIBLIOGRAPHY 


1. M. L. Barry and P. Olofsen, "Advances in Doped Oxides as Dif- 
fusion Sources," Solid State Technology, Vol. 11, No. 10, 

pp. 39-42, October, 1968. 

2. M. L. Barry and P. Olofsen, "Doped Oxides as Diffusion Sources," 

J. Electrochem, Soc. : Solid State Science, Vol. 116, No. 6, 

pp. 854- 860, June, 1969. 

3. R. M. Burger and R. P. Donovan (Editors,), Fundamentals of Silicon 
Integrated Device Technolo gy, Vol. I, Section II, "Diffusion' 1 

by A. M. Smith, Prentice-Hall, Englewood Cliffs, N. J. , 1967. 

4. K. M. Busen, W. A. FitzGibbons, W. K. Tsang, "Ellipsometric 
Investigation of Boron-Rich Layers on Silicon," J. Electrochem. 

Soc.: Solid State Science, Vol. 115, No. 3, pp. 291-294, March, 1968. 

5. T. L. ChuandG. A. Gruber, "Doped Silicon Films, " Electro 
Chemical Technology, Vol. 5, No. 1-2, p. 43-46, Jan-Feb, 1967. 

6. C. S. Fuller and J. A. Ditzenberger, "Diffusion of Donor and 
Acceptor Elements in Silicon, " Vol. 27, No. 5, pp. 544-553, 

May, 1956. 

7. N. Goldsmith, J. Olmstead, and J. Scott, Jr., "Boron Nitride as 

A Diffusion Source for Silicon, " RCA Review, pp. 344-350, June, 1967 

8. A. S. Grove, Physics and Technology of Semiconductor Devices , 
Chapter. 2, (Thermal Oxidation), Chap. 3 (Solid State Diffusion, ) 

John Wiley and Sons, New York, 1967. 

9. A. S. Grove, O. Leistiko, andC. T. Sah, "Redistribution of 
Acceptor and Donor Impurities During Thermal Oxidation of Silicon," 

J. Applied Phys. , Vol, 35, No. 9, pp. 2695-2701, September, 1964. 

10. M. S. R. Heynes, "Boron Diffusion into Silicon Using Diborane, " 

Electrochemical Technology, Vol. 5, No. 1-2, pp. 25-33, January- 
February, 1967. 


40 



11. Shiro Horiuchi and Jiro Yamaguchi, "Diffusion of Boron in 
Silicon Through Oxide Layer," Japanese Jour, of Applied Phys., 

Vol. 1, No, 6, pp. 314-323, December, 1962. 

12. John C. Irwin, "Resistivity of Bulk Silicon and of Diffused Layers 
in Silicon, " Bell Systems Tech. Jour., Vol. XLI, No. 2, pp. 387- 
410, March, 1962. 

13.. T. Kato and Y. Nishi, "Redistribution of Diffused Boron in Silicon 
by Thermal Oxidation," Japanese Journal of Applied Phys., Vol. 3, 
pp. 377-383. 

14. D. P. Kennedy and P. C. Murley, "Impurity Atom Distribution 
from a Two-step Diffusion Process, " Proc. of IEEE, Vol. 52, 

No. 5, pp. 620-621, May, 1964. 

15. James S. Kesperis, "The Effect of Gas Concentration in the 
Diffusion of Silicon from a Phosphine Source, " J. Electrochem. 

Soc.: Solid State Science, Vol. 117, No. 4, pp. 554-557, April, 1970. 

16. A. D. Kurtz and R. Yee, "Diffusion of Boron into Silicon," J. Applied 
Phys., Vol. 31, No. 2, pp. 303-305, Feb., 1960. 

17. K. Lehovec and A. Slobodsky, "Diffusion of Charged Particles Under 
Consideration of a Built-in Field, " Solid State Electronics, Vol. 3, 
1961. 

18. W. R. Runyan, Silicon Semiconductor Technology, McGraw-Hill, 

New York, 1965. . 

19. Guenter H. Schwuttke, Semiconductor Junction Properties as 
Influenced by Crystallographic Imperfections , Report No. AFCRL- 
67-0564, prepared for Air Force Cambridge Research Lab by IBM. 

20. W. M. Whittlemore and G. M. Oleszek, The Effects of Process 
Variables on the Open-Tube Diffusion of Boron Into Silicon From 
Boron Tribromide, IBM Report, TR 01.1111, December, 1968. 


41 



APPENDIX! 


NORMALIZATION AND EXPLICIT DIFFERENCE APPROXIMATION 


The diffusion equation is normalized so that it can be written in 
terms of dimensionless time and space coordinates. The normal- 
ization factors are given in the following two equations; 


Vdt 


( 6 - 1 ) 


T= 


_t_ 

T 


( 6 - 2 ) 


where T is the total diffusion time. 

Writing an equation with this type of normalization facilitates 
the choice of increments to be used in a numerical solution. 

The equation is normalized by using the chain rule for differ- 
entiation. As an example, the time derivative is rewritten as 
follows; 


3N _ _1 3N 

at ~ r 


(6-3) 


The normalized diffusion equation is given in the following 
equation. 


3N 
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(6-4) 



The time and space derivatives are approximated by finite 


differences. 
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( 6 - 6 ) 


The superscripts i+1, i, and i-l indicate that N is evaluated 
at times T+ AT, T, and T-AT respectively. The j subscripts are 
analagous for the spatial coordinate. 

Substituting the finite differences into the diffusion equation 
gives the following difference equation. 
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(6-7) 


Vdtay 


dt 


In equation (6-7) all unsuper scripted variables are assumed to 
be evaluated at time T. 

With the first and last values of N given from the boundary 
conditions, the difference equation can be used iteratively to predict 
the impurity profile at the next time instant if the present profile 


is known. 
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APPENDIX II 


IMPLICIT DIFFERENCE METHOD 


The implicit method is similar to the explicit method except 
that the spatial derivatives are written in terms of the impurity- 

concentration at the next instant of time. 

Ni+I ivi-tT 
BN V 1 - N j-1 


BY 


± 

2 AY 

i+1 ..I 


0 N N l:i - 2 Ni -t-Nj 
9 Y 2 A y 2 
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(7-2) 


Substituting the finite difference derivatives into the diffusion 


equation gives the following difference equation. 
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(7-3) 


In equation (7-3) there are three unknowns, i.e., the three values 
of impurity concentration. If j runs from zero to J, where JAY is 
the largest value of Y to be considered, then there are J-2 equations 
in J-2 unknowns to be solved simultaneously. These equations can 
be written in the form of a matrix equation as is shown in equation 


(7-4). 
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In this equation A, B, and C are the coefficients of N. ., N , 

J‘ L j 

and Nj._^ respectively; and they vary in both time and space. Since 
the coefficient matrix is tridiagonal, equation (7-4) can be easily 
solved by Gaussian elimination for all of the next values of impurity 
concentration. The first and last elements of the column vector on 
the left-hand side of equation (7-4) each contain two terms. The 
boundary conditions determine the second term in both cases; and 
since the boundary conditions are included in the matrix, the 
solution to equation (7-4) will depend on these boundary conditions. 
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APPENDIX III 

IMPLICIT PROGRAM LISTING 

C PROGRAM For CALCULATING impurity profile ^JUNCTION depth# 
c and sheet resistance FOR Two step diffusion in silicon 

DIMENSION N(200> # A (200) »R(200> »E<200> »DA(200) » CON (200) 
REAL NO»LA»Ni3»M 
REAL**) NJ» ANMlf ANPi'AU 

Nl(T)=EXP(21.2S)/( (300. ())**» 1-5) )*l#5El0*KXP(-<6420« 
*-2.O0*T>/T>*T**(i.5> 

RA(N3»N4)=N3/2.+S0RT(N3**2/4.+N4**2) 

C Nb IS THE BACKGROUND DOPING » PER CC 

Nt3=1.0Elb 

C PREDEP TImE(MIN) »Tt:MP(DEG»C. ) rLRlVE-IN TEMPtDEG.C. ) » 

C TIME IN W2(MIN> #02<MIN> »ANS1 STEAM<MIN> 

C PD IS 1.0 FOR PREDEP. U.l)(OK BLANK ) OTHERWISE 

C TOT IS ToiAL DHIVL-IN TIME IN MINUTES 

40 READ (5»10)TI1»TE1» TE2 » T 1 t T2 ♦ T3 ♦ PD » TOT 
1(1 FORMAT (BF 10.3> 

IF (F’O.GT . 0 • ) XK = 0 . 0 
WRITE ( f,» 301 > Til # TE1 »TE2 » T1 1 T2 i T3 
301 FORMAT (' *»6EI0.2) 

IF (Til .LE. 0 .0 ) STOP 

Ul=tn 0 *FXP < -42 . £3/ (TEX + 273.0) ) 
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TI1=T11*60.0 

TOT=TOT*6G.O 

T2=T2+T1 

T3=T3+72 

T1=T1*60.0 

T2=T2*60.C 

T3 = T3*;:0 « G 

C calculation OF OXIDE GROWTH rate parameters 
CO:: 0«i71*EXP(-23.1E3/<TE2+273. ) ) 

CS=2 • 37*cXP(-22.7E3/ ( TE2+273. ) ) 

BS=0.6:>£E-9*EXP(-6 . 2E3/(Tfc' 2+273 . ) > 
fj0=2 • i, •( ~9*EXf 1 ( -14 *4£3/ ( TF.2+273 • ) ) 

T 12- 7 J 

C DIFFUSION CONSTANT CM**2/StCU;oR0N GIVEN) 

D=.4*E:iP(-40.5E3/(TE2+273. ) > 

IF(PD.G1 . .91)TI2=01/D*TI2 
IF (T 12. Lt. 0.0) GO TO 6 
L A— 2 » U * S OR T l D I * 7 II) 

NO— 4 •t.k'C 
0 = NO*L/ 1 .* . bO'i 
ALPHA=!J.'43 

C GAUSSIAN ESTIMATION OF JUNCTION DEPTH 

XJ= ( 4.*D«T0T*AL0G(Q/(ND* <3. 14*D*T0T)** < .5) )))**< .5) 

C CALCULATION OF NORMALIZATION FACTORS 

XO=SGRT (u*TOT) 

DX=3.*XJ/ < 100. *XO) 

L»T=DX**2 
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rt=dt*tot 
RX=OX*XO 
M = IF IX < 1 »/0T ) 

M=M*TI2/T0T 
J=IFIX(3.*XJ/(X0*0X) ) 

X=0. 

1F(F0.LT. 1.0>GO TO 300 

c set initial condition from complementary error function 

oo 1 I=l»j 

N< I )=N0*L'KFC<X/($SRT(01*T1 )*2.C> > 

X=X+UX *Xo 
1 CONTINUE 
60 TO 40 

300 CONTINUE 

L=2 

R=PT/(UX**2) 

DO 3 L~1 rbi 

C DETF.RM3NA 1 ION OF AMR I ENT CONDITION 
1FU.*R r*LT«T3)C0X=CS 
I F (L*RT *LT«T£) COX-CO 
IF (L*RT.uT*T3)B0X=US 
IF (L*RT.L V.T2)00X=ii0 
0 x o = Li 0 X / < O 0 X / C O X ■ + 2 . * X K ) 

IF (L*RT.L7.TI)DXO=0. 

XK=XK+OXO*RT 

J0=J-1 

N(J)= 0.0 
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C SET aOUMuARY CONDITION ACCORDING TO AMBIfc.NI CONDITION 
N ( 1 ) =M ( ?. ) / (RX* ( 3 . 3-AU’HA ) *OXO/D+l . 0 ) 

1 F ( PD • G T • (i • 95 ) N ( 1 ) "NO 
DO 2 I=2*tiO 
GA= • 45*R T /RX«DXO 
ANMl=Ol.LE(N(l-l ) ) 

ANPI=Oi 1+1) ) 

T=TE2 

IF (PD.O i o(j* ) I " f t. .i 

C CALCULATE INTRINSIC CARRIER CCuClNTRATION 

A0=f.;I ( 1+273,0) 

AO=NI ( TC2+273 • 0) 

A0=DbLr.(A0> 

N3=DULf (»4<I> > 

C CALCULATE HOLD CONCENTRATION 

RAN1=RA < AnMI f A 0 ) + 1 , 0 
HAN2=RA C AN PI t AO) +1 . 0 
K AN-RA ( !'.‘3 r AO ) + 1 • 0 
IF(L.I-T.2)RAN=1.0 
IP ( L • L T » 2 ) R A N 1 ~ I • 0 
IP (L »L ■ • ti ) R A!\i2~ 1 • 0 

c calculate terms for matrix Lou at ion 

Ad) =-K/4. * (RAM2-RANJ ) /RANI-GA/2 . 0 + f< 

B< 1 ) =-l • u-2 . *R*K* ( RAN2-2 . *RAN+RAn1 ) /RAN 
L d ) =R/A . 0* ( RAN2-RAN1 ) /RAN2+GA/2 .O+R 
DA(I)=-MD 
? CONTINUE 
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C CALL SUBROUTINE FOR SOLVING MATRIX EQUATION 
0A(2>=DA(;?)-A(2>*Na> 

DA (JO) -f)A ( JO) ~£ < JO) *N< J) 

CALL T!\'IDAG<2» JO»A»BrE #OA»N> 

3 CONTINUE 
GO TON 0 

4 CONTINUE 
6 CONTINUE 

DO 5 1=2 »j 
C=N(I“J ) 

C G<C) FITS IRVINS DATA FOR CONDUCTIVITY 

CON( I-J >=G(C> 

Y=I*RX-2.*RX 

Nl=I-I 

WRITE(0»99)M(I-D*r 
99 FORMAT (• * f 2E20»t>) 

C FIND JUNCTION DEPTH 

ZF (N C I ) .LT.NDJOO TO 30 
5. CONTINUE 

C CALCULATE SHEET RESISTANCE 

30 R-RS (CON* Nl t RX ) 

TJ=T3/fi0.0 

Tl=Tl/60. 

T2=T2/60. 

T3=T3-T2 

T2=T2-T1 

WRITE (6»20)R»Y»Ti»T2rT3*XK» AO 
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20 FORMATt* »»♦ ’/»7£15.3) 
GO TO uO 
END 


real function Gtcx) 

CN= ABS ( Ca ) 
iFtCX.LT.O.O) GO T03 
IF (CN.GT.u.O) A-i « 0 
IF(CN.GT.U*0) B=7.2E-17 
IF ICN.GT.i.bii +16) A=0.b5 
IFCCN.GT.1.5E +lo)B=3.3E-ll 
XF(CN'.GT.i*4E+18) A=0.332 
I F < CN . G T . fc • 4 E-* 1 ti > 0= 1 . 47E- 14 
IF (CN.GT.i •5E‘* 19) A=0.9b6 
IF (CN.GT.l.JjE + 19)b=4*t-l7 
GO TO 6 

iF(CN.CT.u-O) a- 1.0 
IF ( CM ♦ GT • 0 • 0 ) R=2 . E~16 
IF (CN*GT • 3»6E15) A=0«fl37 
I F l C N . G T . 6 • 5E I S ) J=b . 97F> 1 4 
IF(CN.G1 .1 • 0E1 7 ) A=0.54o 
I F { CM . GT . .V • OF 1 7 ) ii-6 • 93E-9 
IF ( CM . GT . 9 . EC 13 ) A=0 . 94 
IF (CN.GT .y-5Elti)U=2.E-lb 
IF (CM » GT * u • E19) A=0.744 
IF ( CM .GT . ti . El9)rl=1.43£-12 
IF ( CM . GT . c . 35E20 ) B= 1 . 04£-t> 
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IF(CN.GT*2.35E20)A=0.45fe 
5 G = B*CCN**A> 

RETURN 

end 


funct ton rs (con # n . nx > 

DIMENSION CON«N) 

SUK-U o !.;• 

X~G • G 

DO 1 I ~ j. # jxj 
B=4.0 

IF ( 1 /. 1 v 2 • t_G • I ) /.)— 2 • C 
IF ( I *T -j . 1 } 'd~l . 0 
IF < I. N) Or- 1.0 

SUM=S!>+C«;N( T.ltP 
1 CONTINUE 

R=f>X/3.*SUM 
RS=3 • /j< 

RETURN 

END 

bUC-ROUT 1 cii_ TR I DAG C JF # L * A » B » C f 0 » V ) 

DIMENSION A (200) »R(20P> «C(2U0) »D <200 ) # V ( 200 ) 
DOUBLE PRECISION bLTA <200 > r GAMMA (200 ) 

ULTA( JD=uUF) 

GAMMA ( JF ) =D ( JF } /BETA ( JF ) 

JFP1=JF+I 
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DO 100 I=jFPirL 

BETA ( I ) -U II ) -A ( I ) *C (1-1 > /BETA ( 1-1 > 

100 GAMMA ( I ? = i 0 • l ) - A < .1 ) *GAKMA < I~1 ) ) /BETA ( I ) 

LAST-L-JF 
V CL) -GAMMA (L) 

DO 200 K-l »LA5T 
I=L-K 

200 v ( i ; -gamma ( i ) ~c c i ) *v u*i) /bet acd 

RETURN 

END 


function ERFC(A) 

0 1 MENS l OH Z < 1 5 ) * n C X 5 ) 

DOUBLE PRECISION X 

HEAL Z /. 09331* .49269' 1 «215for 2.26994 » 3 .6672 * 5.42533 1 
*7 . 5669 » 1 0 , 12022 » 13 . 13026 rib. 6544 r 20 . 77647 » 25 , 62389 » 

$ 31 * 4 0 7b 1 r 36 » 530 66 r 48 « 026 OB/ 

REAL to / .21623. , 34212 ' . 26302 r . 12642 * .40206 r .056 36 1 
••5=. 12124 » . ll 167 * c 64599 f » 22263 » .422 74 r .39218*' .14565 
*r .14865 » .16005/ 

SUM=0 . 

DO 1 1=1 r 15 
X=?UHA 

S U H = S U . + 4 { I ) * L X f > ( - X * * 2 * X ) 

1 CONTINUE 

E RFC- SUM *£XP ( - A ) *2 . 0 / S Qh T (o . 14.159) 

RETURN 


END 
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APPENDIX iv 

EXPLICIT PROGRAM LISTING 


C PROGRAM FOR CALCULATING IMPURITY PROP II F * DUMP T JON DEPTH# 

C AND SHEET RESISTANCE FOP T*’C STEP DIFFUSION IH SILICON 
D I MENS I ON H C 100 ) r NOT ( 1 0 0 ) * CON ( 1 0 0 > >YP l 2(1 0 > i> 7 f 2 [) 0 > 

REAL N»NU« * NO* LA* NO 

C If F'LwT Of iJlA) lb DESIRLU f LAKE TrtF FIRST CARD 0UTPT=2, 

C OTHER C iSL MAKE THE CARP <m JTMTij , a 

outp r~‘! . tj 
OUT PT~, « •;} 

I F ( OUT PT * l_T . 1 • 0 > WRITE <6# 100) 

100 FORMAT (/* • TC1=PR£UEP TEMP « DEG.C.*TEa=ORIVt>IN TEMP. UEG.C.*/ 
•*»RS1=PRF0LP SHEET RESIST ANCE #nS2=DR IVOIN SHEET Pf-.S ISTaHCP. o » 

* ♦ GHNS * / * Xu= JUNCT I ON Dh'Pl Hr Ml. CRONS » XO-OX I « >E Tf i rr r. ; ? c,c; ♦ 

** IN ANGaTR0wSV*TCN2)= TIME In NITROGEN* T(<>2) v T IMF* 

*♦ IN 0XYG£N»/»TIH20) = TIME IN STEAM r ALL IN MINUTES’//# 

*1 OX* ’ TlI * *?X* »TE2**7Xr »KS1 * »7Xr ♦ftS?.* »7Xr fC-Xf ' X 0" f 

*6X * * T ( ,2 ) * * 5X » * T ( 02 ) * » SX * ♦ T ( H20 ) » / ) 

C N13 IS THE BACKGROUND DOPING * PEP CC 
NB=1.0E1d 

C PKEDEP TIME (MINI *TEMP(DEG.C. ) r DRIVE- IN TEMP ( DEO * ( . » ? 

C TIME IN W2U.-IN) *02(MIN) » AND ST r AM (M IN) 

C INITIAL VALUES OF DATA 
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C SET ICON .LT • 2 FOR A BORON NITRIDE PREOEP 

IC0N=2 
T1=0.0 
T2=0.0 
T3=0.0 
IOAT=l 
TU=17. 

TE1=980. 

TE2=1050. 

RS1=72. 

c oata generation program 

40 CONTINUE 

IF(IDAT.GT.6) GO TO 504 

T2=0.0 

T3=0.0 

Tl=Tl+20. 

IOAT=IDAT+l 
GO TO 510 

504 lFdOAT.GT.12) GO TO 506 
U=0.0 
T2=T2+20. 

T3=0.0 
IDAT=IDAT+1 
GO TO 510 

506 IF(IDAT.GT.18) GO TO 508 


T1=0.0 
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T2=0.0 

T3=T3+20. 

IDAT=IDAT+1 
GO TO 510 
508 TE2=TE2+50 • 

IF(TE2.GT.1200.) STOP 
I0AT=1 

GO TO i4 0 
510 CONTINUE 

C DIFFUSION CONSTANT CM**2/StC (HURON GIVEN) 
T-Tt 1+273* 

cond=o. 

GO TO 200 

198 01=2. *(X 
T=TE2+273. 

• C0N0=2. 

GO TO 200 

199 U=OC 

GO TO 13 

200 IF(T-1<423.) 201 » 20 If 202 

201 TACT=ALOG(7.04)*1273.*1423./150. 
DOB=O.GS2-13*EXP < TACT/1273. )/(3*b> 

GO TO 203 

202 TACT=AU06(6.299)*1423**1473./50. 
DOR=0.65E-12*EXF’( TACT/1473. > 

GO TO 203 

203 DC=D0H*EXP (-TACT/T) 
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IF(C0N0.UT.1. ) GO TO 198 
GO TO 199 

13 LA=2.0*SQRT<60.*m*TIl) 

IF ( ICON-21 14 » 15» 15 

14 ALNRS=18. 19-1. 38E-?.*TF 1-6. 24E~2*TI 1+7. E-4*TI1**2 
R5 lrEXP ( ALNKS ) 

QSL=0.n 
G'S 1 = 0.0 

CALL FiNUNO(RS1*LA»NO» 

GO TO 16 

15 CALL FlNONO(RSl#LAfNO> 

QSlsO.O 

C SET QSL=-i.O IP GLASS INTERLAYER IS REMOVCO* 

QSL=3 . bEib-0 . S64«NO*LA 
QSLs-1.0 
18 T2=T2+Tl 
T3=T3+T2 
Tl=Tl*oO»u 
T2=T2*r.0.o 
T3=T3*bO *0 

C OXIDATION RATE STEAMCCS AND BS> AND 02(C0 AND RO) 
CO=0 . l7l*tXP 1-23 . 1E3/ ITE2+273. > ) 
CS=2.37*LXP(-22.7E3/<TE2+273.) ) 
60=2»18E-9*EXP(-i4.4E3/<TE2+273.) ) 

BS=o .6b2E-9*EXP (-6.2E3/ (TE2+273. ) ) 

C TOTAL TlMt FOR GAUSS I AM ESTIMATION 

TI2=T3 
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Q=N0*LA*.564 

ALPHA=C.45 

C GAUSSIAN ESTIMATION OF JUNCTION DEPTH 

C diffusion EONS, integrated OVER REGION OF 2*XJ 

XJ=<4.*D*TI2*AL06(Q/<NB*(3. 1’’-*D^T12 ) ** ( . 6> ) > 
c normalizing factors 

XO=SQRT ( J*TI2) 

QX-2 • ♦XJ/ (200 • *X0) 

DXz2 • */J/ (iOO.+XO) 

DT=DX**2/(4. ) 

RT=DT*TI2 

RX=OX*;<0 

R=DT/lOX**2> 

M=IFIX(1./0T) 

J=IFIX(2.*XJ/(X0*DX> > 

x~o. 

c initialization OF CONCENTRATION 
DO 1 1 = 1 rU 
NU)=N0*EKFC(X/LA) 

X=X+QX*XC 
1 CONTINUE 

CSS=5 . c'Elb* « TE2 ) ** < 1 . S9b) 

GSU=0.S64*NO*LA 

X02=0.0 

WR I T E ( o f 0 0 0 ) X J » D » R X t QS L » QSU t C SS 
300 FORMAT ( / r bE 1 0 « 3 ) 

L=Z 
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21 CONTINUE 
NDT(J)rN(J> 

C BOUNDARY CONDITION AT Sl-SI-02 INTERFACE 
IF(L*RT-T1> 120 » 120 1 122 
122 IF(L*RT-T2> 124»124»126 
C OXIDATION AND DRIVE- IN WITH STEAM 

C IMPURITY FLUX AT SURFACE ASSUMED CONSTANT UNTIL GLASS IS 

C REMOVED OK TRANSFORMED 
126 TIMES=A83(L*RT-T2> 

XK=CS*( ( 1 • + 4 • * (CS**2 1 ** T ImES/BS+4* * i CS**2 ) * ( Xo2**2 ) / (B5**2 1 
*+4.*CS+X02/BS>**( .5>-loO>/(2**CS> 
if(qsl.lt.o.o) go to 142 

IF (XK-B.E-6) 140 r 140 » 142 
140 DXO=C • 

GO TO .1 6 U 

142 DXO =C5/(1«+4.*(CS**2)*C TIMES ) /tiS+4. * ( CS**2 > * C X02**2 ) 

*/ < B S * * <1 ) +4 . *CS*X02/BS ) ** ( . 5 ) 

60 TO 16 Ll 

C OXIDATION AND DRIVE-IN IN DRY OXYGEN 

C IMPURITY FLUX AT SURFACE ASSUMED CONSTANT UNTIL GLASS 

C LAYER IS TRANSFORMED 

124 TlMEO=ABS ( L+RT-T 1 ) 

X02=BC*( < 1 • +4 • * ( C 0 ** 2 ) *T IMEO/BO +8.E-6*Co/80) ** < )/ 

* ( 2 • *C0 ) 

IF {QSL.LT. 0.0) GO TO 146 
IF (X02-0.L“6> 144 r 144 r 146 
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144 DX0=0.0 

GO TO 160 

146 DXO =CO/(l.+4.*<C0**2>*< TIMEO ) /B0+8.E-6*C0/B0> *♦ ( • 5) 
GO TO 160 

C DIFFUSION IN NITROGEN - ACCOUNTS FOR BORON GLASS NOT 

c removed from the surface by the hf etch 

120 lFtQSL.LT. 0.0) GO TO 152 
IF(QSI-QSL) 150 » 150 » 152 
150 F=D*(CSS-im<2> >/RX 
D0SI=F4RT 

qsi=qsi+uqsi 

QXO-0 « 0 
GO TO 162 
152 F=0.0 

DXO=0.0 
GO TO 162 

160 F=(ALPHA~3.33)*0xO*N(i) 

162 GA=0.4b*RT/RX*DXo 
JO=J~l 

C THIS IS AN ARTIFICIAL BOUNDARY CONDITION TO MAKE N(X) 

C FINITE AT 2*XU And EQUAL TO GAUSSIAN VALUE 

IF ( L . GT . 10 ) NOT ( J ) =0/ ( 1 . 79*SGRT ( D*L*RT ) ) *gxp < - < J*RX 
*(4.*D*L*RT) 

C INTEGRATION OF NORMALIZED DIFFUSION EQUATIONS 

DO 2 I -2 » JO 

NOT ( I )=R*M I + D+ (l.-2.*R+GA)*N( I ) ♦ (R-GA > *N < 1-1 ) 

IF (NOT ( I ) ,LT. 1 .OE-15) NOT ( I ) =1 . OE-15 
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2 CONTINUE 

IF ( N tl ) + ( DX*XO*F > /D . LT . 0 . ) F=-NDT { 2 ) * 0 / ( RX > 

NOT < 1 ) -Da*F*xO/D+NDT (2 ) 

DO 3 1 = 1*0 
N 1 1 ) =NDT ( i ) 

3 CONTINUE 
lF(L.6E.Mi GO TO 6 
L=L+1 

GO TO 2l 

tf continue 
6 CONTINUE 
K z 2 

DO 5 1-2*0 

C MAKE NEXT INSTRUCTION CARD r C=N(I-1) FOR P-TYPE 

C MAKE NEXT INSTRUCTION CARD rC=-N(I-l) FOR N-TYpE 

C=N< 1-1) 

C G(C) USES IRWIN’S CURVE FORMULAS 
CONU-l)=o<C> 

Y = l*RX-2 • *RX 

YP(k)=Q*EXP(-(P*kX)** 2/(4.+0*TI2) ) /SORT < 3. 14*D*Tl2> 
YP(K-1)=NCI-1> 

Z<K-1)=I*RX 

P=NI 

YP(l> = YPl<:> 

Z<K)=Z(K-1) 

K=K+2 
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iF(Ntl-l) *LT .NB. AND.N* I > »l.T *N( I-l ) )G0 TO 30 
5 CONTINUE 

C INTEGRATION FOR RSr SHEET RESISTANCE 

30 R S2~ R S ( C Of j f N 1 1 R X ) 

IF ( OUTpT e GT « i • 0 ) CAU. rL0T«Z#YP»K-2*NB*Y> 

XK=ES* ( < 1 .+4.* ( CS**2 ) * <T3«T?! ) /FS+4.* ( CS**2> * ( X02**2> / 

* ( us* *2 > *4 . *CS*X02/RS> **i.5)~1.0)/<2« *Cf > ) 

Tl=Tl/60. 

T2=T2/fcG. 

T3=T3/t>0.u 

T3=T3-T2 

T2=T2-n 

IF ( OUTpT « ol « 1 . 0 ) M l X T E < 6 1 2 0 / ' ' "r Y t T 1 » T 2 e T3 1 Xk » Tt r ?. 

20 format (//f iox» *sheet resistance:: • » Fe • 2 // » i ox * * * iunct ion * 

*• DEPTH iu CMS*»E10o3//f 10Xf *TlMt IN NITROGEN IM • 

** MINUTES - * »F9. 1//- JCX* • f'J.ME IN OXYGEN IN * 

** Ml NOTES s « »F9.1//r *C-Xp»TIME IN STEAM IN MINUTES = « > 
*F9.1//»lOX* f OXIlJ£ THTCM'NgSS IN CM = f *E10*2//f 
*10X» ’DRIVE-IK DIFFUSION TEMP IN DEG* C, = « pF9»l> 

XKA=i « L , 8 *■ X K 
YM=1 *E4»Y 

I F ( OUT j >T . i.T .1.0) WR I TE < 6 r 22 > TL 1 1 ft 2 r RSI t RS2 f YM <• X|<A r T 1 r T2 1 T3 
22 FORMAT UtiOXf 9F1 0 . 3 ) 

400 CONTINUE 
00 TO 40 
END 
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c 


SUBROUTINE FINDNO(RSl»LA#NO> 

FINDNO IS BASED ON THE RELATION SHIP FOR RS FOR AN ERFC 


C PROFILE. 2T WAS CALCULATED FROM A NUMERICAL 
C INTEGRATION# AMD THEM FITTED WITH A PIECEWISE 

C CONTINUOUS FUNCTION. 

real no# la 

GAVE=1.0/(RS1*LA) 

alpa=..:>e-io 
AL pO- • 02oc“l'4 
ALPC=.369t.-lb 
BETA=0.64t> 
b£T0=O.642 
BETC=0.95b 
lF(GAVt-lt>5.) (l #6#6 
4 IF ( GAVE. . LT *40 . ) GO TO 8 
AL 06 C = AL OG < G A V F. / AL PB ) / HE TR 
GO TO 10 

6 alogc=alog<gave/alpo/detc 

GO TO 10 

tt AL0GC=AL0G(GAVE/ALPA)/3ETA 
GO TO 10 
10 NO=EXP(ALOGC) 

RETURN 

END 

REAL FUNCTION ERFC(U) 

IFCU-1.0) 3 » 5 » 5 
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3 N=31 

GO TO 15 

5 IF (U-2.0) 7.9.9 

7 N=41 

GO TO 15 

9 lf(U-2.5) 11.13.13 

11 N=bl 

GO TO 15 

13 IF(EXP<-U**2).LT.1.E-15>ERFC=1.E-1S 
IF (EXP(-U**2> .LT.i.E-15)RETUHN 
ERFC=EXP ( U * «U**2 ) / ( 1 . 772*0 ) 

GO TO 17 
15 SU=0.0 

DX=U/(M-1) 

X=O.U 
1 = 1 

1 8 = 2.0 

IF (1/2* 2.EQ.X)B=4.0 
IF(I.Eg.I) 0=1.0 
IF(I.F.Q.N) 8=1.0 
ER= EXP ( 0 . ) 

SU= SU+ER*B 
X=X +DX 
1 = 1 + 1 

IF(IM-I) 2.1.1 
2 CONTINUE 

5=0X/3,0*SU 
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MFC =l,Q . 

CO^TX/^yjr 

Return 
end 


^• 0 * 5 / 1,772 


teAL ^CflON 6(CX) 

CN= Aflsfcx) 

1F 'CX. t T.(i.oi GO r 0 3 

"‘W.GT.o.o, A=i . 0 

if ‘cn.gt. 0 . 0 , e=v . aE . 17 

,F,C "* sr “-* ♦*» *-o. 6s 
^CN.GT.i.* *i 6)BS3 , x _ u 

^r.z.uEne, A=0 . fl32 

W ‘CN.ST. i! ., E+1<|JB81 

, c N . 0 T . l . 5 t + 19) u - 4>t _ 17 
GO To s 

s IF( CN.c.r. 0 .o, A=1 . 0 
"■ | CN '«. 0 . 0 ) 6: a . E ., 6 
IF( W.GT. 3 .SU5)Aro. n 3 7 
,F( CN.GT 0 .5 E15)a=6>97l: _ ii( 

IF(c n.gt. 1 .o E17) ar0tS43 

1 F< CN-sr.i.o E i 7)B=6 . 93£ _ 9 

1F <CN.GT.9.e, Ele) A=0#9 , t 

IF ‘ CN ‘ &T - 9 - 5F l«'Srs. E . J6 

1 f <cn.gt. 6 . e19J 

ei 9>B=1.43 E _ 12 

IV -12 



1F(CN.GT.2«35E20) A=0.456 
J.F (CNsGT •2*35^20) B=l«04E-6 
5 6 = B*(CN**A) 

RETURN 

END 

FUNCTION kS(CON#n#DX> 

DIMENSION CON IN) 

SUM-0 «0 
X=0.0 

DO 1 I=1»N 

B~2 « 0 

IF ( 1/2*2 «t.Q . I ) 8=4*0 
IF (!•£«. 1)13=1*0 
|F(l.EQ.N)Bsl*0 
■ SUM=SU«+C0 N( i>*b 

i continue 

R=DX,/3.*SUM 

RSsl./R 

RETURN 

END 

SUBROUTINE PLOT(XlrYl*N*NB*XJ) 

REAL NB 

DIMENSION XI (250) rYl (250) » ARRAY (45* 101 ) 
DATA AST * D ASH * B AR * BL AUK / 1 M* * 1H- • 1 H 3 • 1 H / 

DO 19 L.-I * 45 


IV- 13 



DO 19 M-l #101 
ARRAY (L»M)= BLANK 

19 CONTINUE 

DO 20 M=l»101 
ARRAY <1»M)=DASH 
ARRAY (45»M) = DASH 

20 CONTINUE 

DO 22 L=JL#45 
ARRAY U. t U=BAR 
ARRAY (L»10l)=liAR 
22 CONTINUE 

CMAX=AES(Yl(l>-NB) 

CMlN=Ai:SCYl(l)-Nrf) 

XMAX=X1(N) 

DO 4 1=2 rN 
Z1=ABS( Yl(I)-NB) 
IF(ZI.OT.CMAX) CMAX=Zl 
IF(Zl.LT.CMIN) GO TO 12 
GO TO 14 
12 CMIN=21 
GO TO 14 

14 IF ( Z1 *i-T • I » 0 ) Zl = 1.0 
22=AL0Gl0(Zl) 
L=(85.~4.*Z2) 

M=(101. *lXl(I) >/XMAX+.5> 
lF(L.LT.l) L=1 
IF(L.GT.4B) L=4t> 
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IF(M.LT.I) m=i 
IF(M.GT. lOl) m=ioi 
ARRAY (L*M)~AST 
4 CONTINUE 

XOM=l.E4*XJ 

WRITE (6» 30 ) CMAXrCMIN# X*JM 

30 FORMAT (lHlr lOXf »CMAX= * »E6.3»5X» «CMIN= •»E8*3»5X* 
*»XJ (MICRONS) =♦ *F6.3/) 

L=1 

WRITE (6* 80 ) (ARRAY (L#M) • M=1»I0I) 

80 FORMAT (* • » 10H — 10XA1 ) 

L-^ 

WHITE (6»82> (ARRAY (L#M) *M=1»101) 

82 FORMAT (« ♦ » 1H3» * C(X) »»101A)) 

L=3 

WRITE ( 6 » 84 ) (ARRAY (L#M) »M=1»101) 

84 FORMAT ( * « rlHUr9Xr 10J.AI ) 

L-4 

WRITE (6* 86) (ARRAY (L#M) rM=l* 101> 

86 FORMAT (’ » » 111 3f * PER CC '»10J.A1) 

IT=1 

DO 32 L=5,4S 
EX=10.**( (L-l)/4) 

CX = I,E£l/t.X 

1F( (L-D .hG. (4*IT> ) GO TO 40 
WRITE (6 r 88 ) ( ARR A Y ( L f M ) » M= 1 r 1 0 1 ) 

88 FORMAT(* « » J.H3f9X»10lAl) 
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GO TO 32 

40 VJRITE<6#90) CX* ( ARRAY IL»M) »M“1» 101) 

90 FORMAT (• »»2Hi2 »£6.1#2H— »101A1> 

IT~lT+i 
32 CONTINUE 

WRITE<6»94> 

94 FORMAT (• * t 1H3»9X#IH3»9X# IH3#9X* IH 3#9X* lH3r9X» IM3»9X> 1H1*9X» 
*1H J»9X»lHJ»9Xf IH3#9X»2HJ»9X> 

XAA=2.E3*XMAX 

XAB=4.E3*XMAX 

XAC=fe*£3«XMAX 

XAD=8.C3*XMAX 

XAE=1.£4*X«AX 

WHITE ’ f : » 96 ) XAA t XAB rXAC » X AD» XAE 
96 FORMAT (/ » 10X» 1H0» l7X#Fb.2» 3 5Xr F*>»2» li>X#F5.2» 15X»F5.2> 
■*15X*F5.2J 
WHITE < f>r 96 > 

96 FORMAT ( /»45X* »X -MICRONS* ) 

RETURN 

END 
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APPENDIX V 


JUNCTION DEPTH, SHEET RESISTANCE, 
and OXIDE THICKNESS CURVES 

The curves given in this appendix are for drive-in diffusions in 
nitrogen, oxygen, and steam at drive-in temperatures of 1050°C, 

1100° C, 1150° C, and 1200° C. The predeposition diffusions are 
assumed to be the 72 ohm or 27 ohm diffusions carried out at 980° C 
that are done at MSFC. The 72 ohm predep is for 17 minutes in 
diborane followed by 7 minutes in O 2 . The 27 ohm predep is for 60 
minutes in diborane followed by 6 minutes in N 2 . All the boron glass 
is assumed to be removed for the 72 ohm diffusion, but the boron- 
rich interlayer is assumed to remain for the 27 ohm diffusion. The 
explicit integration program listed in Appendix IV was used to generate 
the data. Experimental data is given on graph V-5 which has some 
scatter in the resistance after predeposition. All experimental data 
points are for 1150° C diffusion and are shown as circles with data 
spread indicated by a vertical bar. 
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